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ABSTRACT 

In this work, we study a class of early dark energy (EDE) models, in which, unlike in standard dark energy models, a 
substantial amount of dark energy exists in the matter-dominated era. We self-consistently include dark energy perturbations, 
and constrain these models using current observations. We consider EDE models in which the dark energy equation of state 
is at least w m >, -0.1 at early times, which could lead to a early dark energy density of up to Q.de(zcmb) = O.03f2 m (zowfl). 
Our analysis shows that, marginalizing over the non-dark energy parameters such as Q m ,Ho,n s , current CMB observations 
alone can constrain the scale factor of transition from early dark energy to late time dark energy to a, 0.44 and width 
of transition to A, £ 0.37. The equation of state at present is somewhat weakly constrained to wo ^ - 0.6, if we al- 
low Ho < 60 km/s/Mpc. Taken together with other observations, such as supernovae, HST, and SDSS LRGs, wo is con- 
strained much more tightly to wo <, -0.9, while redshift of transition and width of transition are also tightly constrained to 
a, 0.19, A, 0.21. The evolution of the equation of state for EDE models is thus tightly constrained to ACDM-like be- 
haviour at low redshifts. Incorrectly assuming dark energy perturbations to be negligible leads to different constraints on the 
equation of state parameters- wo <, -0.8, a, <, 0.33, A, 0.31, thus highlighting the necessity of self-consistently including 
dark energy perturbations in the analysis. If we allow the spatial curvature to be a free parameter, then the constraints are 
relaxed to w < -0.77, a, ,< 0.35, A, ,< 0.35 with -0.014 < fl K < 0.031 for CMB+other observations. For perturbed EDE 
models, the 2a lower limit on as (os > 0.59) is much lower than that in ACDM (erg > 0.72), thus raising the interesting pos- 
sibility of discriminating EDE from ACDM using future observations such as halo mass functions or the Sunyaev-Zeldovich 
power spectrum. 



1. INTRODUCTION 

Over the last decade, the unexpected faintness of dis- 
tant Type la supernovae have shown th at the expansion 
of the universe is ac ce lerating at presen t (iRiess etal. 1998 : 
Perlmutter et al. 19991 iTonrv et al. 20031: iRiess et al. 2005 ; 
Astier et al. 20051: IRiess et al. 20071 IWood-Vasev et al. 2007 : 
Kowalski et al . 2008llHicken et al. 2009llLampeitl et al. 20091 . 
This remarkable discovery points to the existence of dark en- 
ergy (DE), a negative pressure energy component which domi- 
nates the energy content of the universe at present. Other, com- 
plementary, probes such as the Cosmic Microwave background 
(CMB) and various large scale structure surveys have also con- 
firmed the existence of this mysterious component of energ y 
dEisenstein et al. 20051 iKomatsu et al. 20091 iReid et al. 2 009). 
Several theories have been propounded to explain this phe- 
nomenon, the simplest of which is the cosmological constant 
A, with a constant energy density and a constant equation of 
state w = — 1 . The cosmologic al constant is fit well by the cur- 
rent data dHicken et al. 20091) . however, there are no strong 
constraints on the time evolution of dark energy at present. 
Thus, evolving models of dark energy remain viable as alterna- 
tive candidates for dark energy. Many non-cosmological con- 
stant phenomenological explanati ons for cosmic acceleration 
have been sug gested (see reviews ISahni & Starobinskv 2000; 
Carroll 200U iPeebles & Ratra 20031 ICopeland et al. 2 006; 
Noiiri & Odinstov 20071 ~ ISahni & Starobinskv 2 006; 

Frieman et al. 20081 and references therein). These are based 
either on the introduction of new physical fields (quintessence 
models, Chaplygin gas, etc.), or on modifying the laws of grav- 
ity and therefore the geometry of the universe (scalar-tensor 
gravity, f{R) gravity, higher dimensional 'Braneworld' models 
etc.). As of now, there is no consensus on the true nature of 
dark energy. 

An interesting class of models which have been suggested 
in the literature are early dark energy models, a class of 
dark energy in which the early universe contained a substan- 



tial amount o f dark energy. These models were studied the- 
oretically in (iDodelsonefa/. 20001 ISkordis & Albrecht 20021; 
iDoran & Robbers 20061) and references therein, and have 
been analyzed with respect to observa t ions extensively in 
recent times in dL inder & Robbers 2008; Fran cis et al. 2008; 
iGrossi & Springel 20091 iFedeli et al. 20091 IXia & Vie!"2 009). 
For now, there are no strong observational constraints on the 
EDE models, and it is especially difficult to discriminate EDE 
models which have w = -1 at present from the ACDM model of 
dark energy. 

In this work we use a parameterization of the equation of 
state of dark energy to study and constrain EDE models using 
the currently available data. We attempt to see if bounds can 
be put on the transition from early dark energy to the present 
day dark energy content of the universe. Section|2]explains the 
methods and data used for this analysis, section [3] shows the 
results, and in section|4]we conclude. 



2. METHODOLOGY 

Dark energy perturbations for dynamic dark en- 
ergy models have been studied in a number of works, 
usually under the formalism of a minimally coupled 
scalar fi el d (See \Ma etal. 1999): |H wan g & Noh 20011 
Hu 20021 iMalquarti & Liddle 2002t IWeller & Lewis 2003 
Bean & Pore 20041 iDutta & Maor 20071 iMota et al. 20071 
iNovosvadlvi & Sergiienko 20081 iJassal 20091 and references 
therein). For practical purposes, e.g., analyzing dynamic dark 
energy models in light of data, it has sometimes been the prac- 
tice to consider dark energy perturbations as negligible, and 
sometimes not. As shown in (Par k et al. 20091) . not taking into 
account the dark energy perturbations correctly can lead to 
erroneous, gauge-dependent results. In our analysis, we self- 
consistently include the dark energy perturbations for the EDE 
models, and also show how the results differ if these perturba- 
tions are not included. 
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2.1. Dark Energy Perturbations 

A homogeneous and isotropic large scale universe can be de- 
scribed by the Friedman-Lemaitre-Robertson- Walker (FLRW) 
metric 

ds 2 = a 2 (rj) [drf- + 8 a pdx a dx fi ] , (1) 

where r\ is the conformal time, dx is the length element, and 
a{rj) is the scale factor. The speed of light c is set to unity, so 
that the time variable has dimensions of length. 
First order perturbations take the form 

ds 2 = a 2 (r)) [(l + 2^(x,7 1 ))dT 1 2 -(l+2^(x,r ] ))S a pdx a dx^] , 

(2) 

where $, are the Bardeen potentials. If proper isotropy of the 
medium is zero, then $ = —iff. 

We adopt two equivalent approaches to account for the dark 
energy perturbations, the first consists of considering the dark 
energy component as an additional fluid, while in the second 
dark energy is defined as a minimally coupled scalar field. Both 
approaches lead to the same result within the framework con- 
sidered, and each has its usefulness in analyzing the results. 

2.1.1. Dark Energy as a Fluid 

In this section we follow the treatment of 
(IWeller & Lewis 20031) . Along with the matter and radiation 
components, we consider dark energy to be an additional fluid 
component, so that the dark energy perturbations are character- 
ized by an equation of state, an adiabatic sound speed and an 
intrinsic entropy perturbation- 
Pde 
Pde 
Pde 
Pde 
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Pde wde 

Defining the frame invariant quantity c 2 sj (the fluid sound speed 
in the frame comoving with the fluid), the continuity and Euler 
equations giving the evolution of the density contrast and ve- 
locity of a fluid with equation of state w, = Pi/pi, and adiabatic 
speed of sound c 2 ai = Pi/ pt, may be written as (prime denotes 
derivative with respect to rf) 

5\ = -3U(c 2 j - wdSt - 9U\c 2 si - c 2 a d( 1 + Wi ) y 



-(1 + Wi )k Vi - 3(1 +wd& 

1 kC 2 :Si 
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where A is the acceleration (A = in the synchronous gauge, 
A = — in the Newtonian gauge), and "H = a' /a = aH is the con- 
formal Hubble parameter. For the matter component, w m = c 2 = 
c 2 = 0. For the dark energy component, a fluid with varying 
wde > -1 has c 2 aDE = w DE -[dw DE /d{\n a)]/3(l + w DE ). For 
scalar field like dark energy models, c 2 sDE = 1. For a more gen- 
eral class of models, such as k-essence, c 2 DE could be vari- 
able as well. To reduce the number of parameters, we con- 
sider c 2 DE = 1 in our analysis, which would still allow us to 
study a wide range of dark energy models. Thus for an uni- 
verse containing matter (CDM+baryons) and dark energy, a set 
of four perturbation equations may be defined for the gauge- 
independent variables 6„,,v m ,5DE,VDE and solved using adia- 
batic initial conditions. 



2.1.2. Dark Energy Perturbations in Scalar Field Formalism 

An equivalent and convenient approach for studying the dark 
energy perturbations is to regard the dark energy component as 
a a minimally coupled scalar field Q with self-interaction po- 
tential V (Q). The field dynamics are given by 

Q" + 2UQ' + a 2 ^- 2 =0, (8) 
dQ 1 

and the perturbations of the scalar field evolve through the per- 
turbed Klein-Gordon equation 



MJ" + 2UMJ'+ [ k2 + a2 ^f) 5Q = 4Q'*'-2a 2 ^ . 



The metric perturbations evolve as 



+ 3H^' + 8nGa 2 V * = 4ttG ( Q'SQ' - ^Jq^Q 



,dV 



(9) 



(10) 



The matter density contrast may be obtained from the above 
equations to be- 

1 



S,„ = ■ 



[3H^' + {STrGa 2 (p m + V) + k 2 }^ 



4irGp 

+4kG [Q'SQ' + a^SQ 



(11) 



The fluid parameters for the dark energy component (as defined 
in section [2". 1.1 i are related to the scalar field variables by 



w D e = 



Q' 2 -2a 2 V(Q) 
Q' 2 + 2a 2 V(Q) 
2 _ 2 dV a 2 
" a - DE ~ + 3dQQm 
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- [Sde - 3H( 1 + wde)vde] 



and the gauge-independent perturbation variables by 

Sde 



-^—lQ'6Q' + a 2 ^-6Q-Q' 2 *] 
a L p D E dQ 
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2.1.3. Imprint of Dark Energy on Observables 

The two basic dark energy dependent observables are dis- 
tance and growth rate. Distance measures are based on standard 
candles, rulers, or number densities as a function of redshift; 
growth rate measures are based on density perturbations in lin- 
ear theory. All distance measures are ultimately based on the 
comoving distance to redshift z 



z dz 



da 



/ H(z) J a aH(a) ' 

e.g., the SNe Type I a observations measure the magnitude of 
distant SNe, given by m B (z) = 51og 10 [(l +z)r(z)] + M. The ef- 
fect of dark energy for distance measures is through the back- 
ground expansion of the universe, i.e., from the Hubble param- 
eter H(z) = H{a)/a. For CMB data, this comes in through the 
angular diameter distance and the sound horizon 



Da(o) = a 



s(a) ■■ 



da 



j a 2 H(a) 7i aH(a) 
a c s (a) da f c s (a) da 



" da 



o a 2 H(a) J aH(d) 



(18) 
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The density perturbations are affected by the presence of 
dark energy firstly through the Hubble parameter, and secondly 
through the linear perturbation of dark energy, as in eqs (|6j |7), 
or eqs ITOl ITTb . For the CMB power spectrum, the effect of 
these is felt most strongly in the ISW effect at low Z, as well as 
in a shift of the peak positions. 

The low / observations can be understood as follows. The 
behaviour of the temperature anisotropy power spectrum in the 
CMB is given by the covariance of the temperature fluctuation 
expanded in spherical harmonics 



Q=4k I —V x \Ai(_k,7]o) 



(20) 



where V x is the initial power spectrum, r\o is the conformal time 
today, and A[{k,r)o) is the transfer function at each I. 
On large scales the transfer functions are of the form 

A,(k, Vl) ) = A}- ss (k) + Af w (k), (21) 

where Aj jSS (A;) are the contributions from the last scattering 
surface from the ordinary Sachs-Wolfe effect and temperature 
anisotropy, and A{ sw (£) is the contribution due to the change 
in the potential along the line of sight and is called the inte- 
grated Sachs-Wolfe (ISW) effect. The ISW contribution can be 
written as 



Af w (£) = 2 



(22) 



where r{rf) is the optical depth due to scattering of the photons 
along the line of sight, and j/{x) are the spherical Bessel func- 
tions. 

The frame-invariant potential <j>, defined in terms of the Weyl 
tensor, is equivalent to the Bardeen potential in the absence of 
anisotropic stress and given by the Poisson equation 



-4nGa 2 Sp , 



(23) 



while its derivative in a matter plus dark energy universe, which 
is the source term for the ISW contribution, is given by 



k 2 (f>' = -4ttG— [a 2 (5p~^+Sp DE )] 



(24) 



From the above equations, it is clear that the magnitude of the 
ISW contribution is dependent on the late time evolution of the 
total density perturbations, therefore on the dark energy per- 
turbations. It should be noted however, that these are not in- 
dependent of other cosmological parameters, and the effect of 
dark energy could be masked due to the degeneracy of the dark 
energy parameters with other parameters such as Ho and the 
curvature of the universe. 

2.2. Parameterization of Equation of State of Dark Energy 

To study EDE models under this formalism, we consider a w- 
parameterization wh ich may represent a la rge class of varying 
dark energy models dCorasaniti et al. 20031) 



w(a) = w + (w m -w )- 



l +e a,/A, l_ e (a-l)/A, 



(25) 



\ +e (a-a,)/A, 1_ e l/A, ' 

where wo is the equation of state of dark energy today, w m is the 
equation of state in the matter dominated era, a, is the scale fac- 
tor at which the transition between wq and w m takes place, and 
A, is the width of the transition. If w m is allowed to be a free 
parameter this parameterization can encompass a large class of 
models, including ACDM and w = constant models. Models 
with constant or slowly varying w ~ — 1 would be consistent 
with current observations, however these are not EDE models, 



as they have negligible amounts of dark energy at early times. 
For such models, there would be very poor constraints on the 
transition parameters, since no significant transition takes place 
between early time and late time dark energy. Allowing these 
models in the analysis would therefore cause the constraints 
on a,,A t to weaken. Leaving the amount of early dark en- 
ergy free would be interesting when comparing EDE models 
with ACDM and other dark energy models. Such comparisons 
have previously shown that while it is possible to put an up- 
per limit on the amount of early dark energy, it is not possi- 
ble to put strong constraints on the evolution of dark energy if 
all the d ifferent dark energy mo d els are considere d. Previous 
studies dDoran & Robbers 20061: IXia & Viel 20091) have con- 
strained early time dark ener gy density to ~ 3% of the matter 
density, however, as seen in dXia & Viel 20 09). the evolution 
of dark energy is weakly constrained. In this work, we study 
the EDE models exclusively, to put constraints on the transition 
from early to late time dark energy. If we are able to constrain 
the minimum redshift (or maximum scale factor) at which such 
a transition occurs, we would know that any signature for EDE 
would be found only in observations beyond that redshift. This 
would also put a constraint on the evolution of dark energy at 
low redshifts. For studying EDE models with this parameteri- 
zation, we therefore choose w m > —0.1, to ensure the presence 
of adequate amounts of dark energy at early times, so that we 
may put constraints on the transition from early to late-time 
dark energy for these models. 

2.3. Observations 



We use the latest version of COSMOMC dLewis & Bridle 20021) 
for our analysis, modifying the CAMB module, as well 
as the various modules pertaining to large scale structure 
and supernova observations in COSMOMC, using the equa- 
tions defined in section 12.11 For the analysis using onl y 
CMB data, we use the 5 yr WMAP dKomatsu et aL20 Q9). 
CBI dRedhead et al 2004. VSA dPickinson et al. 20 04). 
BOOMERANG dPiacentini et al. 20061) and ACBAR 
dReichardt et al. 2009b datasets. In addition to the CMB data, 
we use other observations as well. For supernov ae, we use the 
Constitution dataset (SALT) dHicken et al. 20091) . This dataset 
comprises of 397 Type la SNe, of which about 200 are at red- 
shifts z < 0.1, and the remaining are distributed between z = 0. 1 
and z = 1 .7. We also use the lat est SPSS data re lease (DR7) lu- 
minous red galaxy (LRG) data dReid et al. 20 09). and the recent 
value of the Hubble constant from the SHOES (Supernovae 
and Ho for the Eq uation of State) program, Hq = 74.2 ±3.6 
km/s/Mpc (la) dRiess et al. 2009b . w hich updates the valu e 
obtained from the Hubble Key Project (Freedm an et al. 200 lb . 
We incorporate a top-hat prior on the age of the Universe, 
10 Gyr < to < 20 Gyr. The addition of these other observa- 
tions allows us to constrain parameters such as Ho which might 
otherwise be degenerate with the dark energy parameters of 
interest to us. 

3. RESULTS 

We first study the effect of the different dark energy compo- 
nents on the observations. To this purpose we choose two dark 
energy models- (i) a dark energy model with constant equation 
of state wde = _ 0.9 (ii) an EPE model with wq = —1.0, w m = 
-0.1, a, = 0.3, A, = 0.2. We compare the behaviour of these two 
models of dark energy with that of a ACPM (wde = _ 1) model. 
All three models have identical values for the non-dark en- 
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ergy cosmological parameters (e.g., Qo m ,Ho). The first model 
is chosen for comparing the behaviour of non-perturbative and 
perturbative dark energy for a dark energy model close to the 
cosmological constant in behaviour, while the second is cho- 
sen for specifically studying how early dark energy affects the 
results. 

3.1. Effect of Dark Energy Perturbations on Observable 
Quantities 

Following dWeller & Lewis 2003) we look at the effect on 
observations using the dark energy as a fluid framework. We 
first study the effect of non-perturbative dark energy on the ob- 
servations. In an universe containing matter and a smooth dark 
energy component, the matter perturbations may be calculated 
from Eqs (191 HOI 1TB. using SQ = SQ' = 0, to be- 

S£-HS' m -4nGp m 5 m = 0. (26) 
From this equation, we see that the dark energy component ap- 
pears only in the second term which is effectively a damping 
term, therefore a non-negligible amount of smooth dark en- 
ergy would suppress the clustering of matter at large scales. 
Thus the only effect of dark energy for a smooth dark en- 
ergy model arises through the dark energy density, for both ge- 
ometric (e.g., Type la SNe) and perturbative (e.g., CMB, mat- 
ter power spectrum) data. For matter-dominated regime, the 
above equation would result in S m oc a. For the DE model with 
a constant Wde > — 1, the transition between matter and dark en- 
ergy happens earlier than for wde — — 1. and more slowly, thus 
constant wqe > -1 models are expected to have a smaller con- 
tribution to the IS W effect than ACDM. In fig [TJ (a), we show 
the expansion history of the DE model considered as well as 
that of ACDM. We see that the dark energy density equals mat- 
ter density earlier in the dark energy model, and we expect this 
to have a noticeable effect in the scalar C/'s for CMB data. For 
the EDE model, things are slightly different, as seen in figureQ] 
(d). Since the value of the equation of state today is Wq = — 1, the 
transition from matter to dark energy occurs at nearly the same 
time as on ACDM. Also, because initially the dark energy den- 
sity is higher in this model, this transition is flatter. Therefore 
we may expect that these models would have a larger contribu- 
tion to the ISW effect. The effect on the matter perturbations is 
a mild suppression for both cases as expected from eq d26l i. seen 
in figure[TJ(c), (f). 

Setting the dark energy perturbations to zero artificially is 
however not consistent with the general relativity framework 
except in the case of a cosmological constant, wde = — 1. We 
therefore now add the dark energy perturbations to the cal- 
culation. We consider the gauge comoving with dark matter, 
in which the acceleration is zero. If Sde is initially zero, we 
see from eq. © that it is sourced by the other perturbations 
if wde y _ l via the the source term 3(1 + wde)^' ■ An over 
density causes a decrease in the local expansion rate so that 
\&' < 0. In this case a fluid starts to fall into overdensities if 
wde > -1. In the subsequent evolution of dark energy pertur- 
bations, if c 2 sDE = 1, then the source term for the velocity vde 
is positive, thus causing the velocities to be anti-damped. For 
the density contrast Sde, when k <<H, the term (1 + wde^vde 
can be neglected and the velocity and wavenumber enter only 
via the combination (1 + Wde)vde/K which is small. Thus 
the evolution of Sde is almost ^-independent at large scales, 
and the two remaining source terms -3H{c 2 D e ~ w de)Sde and 
-3(1 +wde)^' are of opposite signs with Sde > and < 
initially. Therefore dark energy perturbations change sign at 



very early times and start decreasing, having the opposite sign 
to that of the matter perturbations, which source Sde through 
the now increasing For the DE model with constant wde, 
we see this effect in figure[TJ(b). Thus Sde and S,„ have opposite 
signs, and at late times when the dark energy becomes a signif- 
icant fraction of the energy density, the total density perturba- 
tions are smaller than those without dark energy perturbations. 
So there is a larger overall change in the potential <// in eq ( |2TT >. 
and the ISW contribution is increased. Since the total decrease 
in Sde is small as wde is close to the ACDM value of — 1, the 
matter perturbations do not change significantly, as seen in figOJ 
(c). For a DE model with wde >> — 1, the effect on S m would 
be stronger. 

For EDE, varying wde provides a further effect. Initially 
when w m =—0.1, the source terms approximate to —3^f' and 
-3HSde, and since iff' is significantly larger than Sde, it is 
the primary source term in eq (|6|. Therefore, Sde decrease 
rapidly, more than it would for wde ~ — 1. When the dark en- 
ergy equation of state transitions from w m = -0. 1 to wq = -1 .0 at 
a, = 0.2, the source term -3'H(c 2 de -Wde)Sde becomes larger 
and therefore the decreasing Sde starts to increase, though not 
fast enough to change signs again, as seen in figure[TJ(e). Sde is 
therefore still of opposite sign to 8 m , but less negative than for 
a wde = constant ~ — 1 case. Thus the ISW contribution is de- 
creased from what it would be in the no perturbation case, but 
still is larger than that for the ACDM model, while the matter 
perturbations at low k, which source the dark energy perturba- 
tions through 3(1 +Wde)^', become smaller at late times as wde 
becomes more negative. Thus matter perturbations at low k for 
EDE models are strongly suppressed at late times as compared 
to ACDM, or the no perturbation case (figQ](f)). The change in 
potential <f>' in eq ( |2TT ) is therefore enhanced. So effectively, we 
expect a strong enhancement off the transfer function and there- 
fore the matter power spectrum at large scales (low k). Thus the 
matter power spectrum at late times, when normalized at low k, 
would show a strong suppression on the small scales (i.e. at 
high k), and this suppression is effected due to the variation of 
the dark energy equation of state. 

The effect of dark energy perturbations can be understood 
also from the scalar field formalism. From eq (0, the scalar 
field Q can be viewed as a fluid with comoving Jeans mode 
given by the curvature of the potential, i.e. the mass of the 
field, kj = a\J d 2 V / dQ 2 . Therefore scales which corresponds 
to modes k < kj will collapse under gravitational instability, 
while modes k > kj will undergo a series of damped oscilla- 
tions due to pressure waves in the quintessence fluid. This has 
two major effects. Firstly, the large scale clustering of dark en- 
ergy enhances the amplitude of the ISW effect in CMB at low 
/. Secondly, as a consequence of the homogeneity of of the 
dark energy component on small scales and the fact that the 
growth of the linear matter perturbations is suppressed due to 
the lower values of ftom, the linear matter power spectrum at 
small scales will have an amplitude which is smaller than in 
ACDM. We thus expect that on the very large scales (k < kj ) 
the dark energy clustering enhances the matter power spectrum 
compared to the unclustered case, while on small scales (k> kj 
) the opposite occurs. If we CMB normalize the matter power 
spectrum (i.e. normalize it at large scales), the small scale mat- 
ter power spectrum will show a stronger suppression of power 
than in the no perturbation case, thus giving a smaller value of 
a% at present. 

Fig|2](a), (b) show the CMB C/'s and the matter power spec- 
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trum at z = normalized to CMB for the DE model. As ex- 
pected from the arguments in the previous paragraphs, we see 
that there is a slight shift in the CMB peak position as well 
as enhanced power at low / for the DE model as compared 
to ACDM. The main effect is at low I, a region which is cos- 
mic variance limited, therefore difficult to rule out observation- 
ally. For the matter power spectrum, as expected, there is a 
small suppression of power at high k (since the normalization 
is done at low k). The value of a% in the no perturbation case is 
a% = 0.79, while that in the perturbed case is <j% = 0.80, and that 
for ACDM is a% = 0.82. Neither the effect on CMB nor that on 
the matter power spectrum is in itself good enough to rule out 
the DE model, even for the case where DE perturbations have 
been accounted for. For the EDE model, as seen in fig [2] (c) 
(d), the non-perturbative case shows effect mostly in the low I 
regime through the ISW effect, which is cosmic variance lim- 
ited. The results for the matter power spectrum today also show 
a very slight difference from the cosmological constant. These 
results appear to suggest that just the non-perturbative effects of 
dark energy are not sufficient to discriminate this EDE model 
from ACDM, especially if we factor in degeneracies with other 
cosmological parameters, such as Hq. When we consider the 
perturbative case, the ISW effect is actually muted, however, 
there is a slightly larger shift in the CMB peak position, (see 
inset of fig[TJ(c)) which is a tightly constrained observable. The 
matter power spectrum at present shows a stronger suppression 
at small scales which leads to a much smaller value of erg = 0.69 
(as compared to the non-perturbative case, where a% = 0.81, 
which is close to the ACDM value). Thus, although the back- 
ground expansion of this model is very similar to ACDM at late 
times, its early time behaviour leaves signatures for discrimi- 
nating it from the ACDM model provided the dark energy per- 
turbations are accounted for properly. The effect of adding 
the dark energy perturbations is seen in fig [2] (e), (f) for both 
DE and EDE models. In obtaining the scalar C/s, for the DE 
model, there is a fairly large difference at low /, while at high I 
the perturbed and non-perturbed models behave similarly. For 
the EDE model, there is a large difference at low /, and also 
a significant difference at the higher Is. For the matter power 
spectrum today, the EDE model shows a larger difference in in 
the perturbed and non-perturbed case. Thus, a model close to 
ACDM today as also in the past (as in the DE model chosen) 
would be difficult to discriminate from ACDM from current ob- 
servations, but a model with a different expansion history in the 
past, even if it is very similar to ACDM today (such as the EDE 
model), could be discriminated using the perturbative observa- 
tions such as CMB and the matter power spectrum provided the 
dark energy perturbations are not neglected. Thes e results are 
commensurate with those found in dMa et al. 1999b where con- 
stant equation of state models o f dark energy were considered, 
and those in ( Alimi et a l. 20091) . where quintessence models of 
dark energy were studied. 

We note here that, since in addition to the ISW effect, dark 
energy also makes itself felt in a shift of the CMB first peak 
position, we expect that the dark energy parameters may be de- 
generate with VI Ji 2 if the flatness condition is removed in the 
analysis. We study the effect of curvature on the scalar C/'s in 
figurefj] A non-flat ACDM model will differ from a flat ACDM 
model with all other parameters identical mainly in a shift of 
the peak positions. Figure[3]shows this shift for a ACDM model 
with Q K =0.06. An EDE model with wq =-0.65, w,„ =-0.1, a, = 
0.2, A, = 0.1, and a curvature Sl K = 0.06 is also shown. For the 



EDE model, the dark energy component compensates for the 
curvature of the universe, thus the peak position is the same as 
for the flat ACDM model. However, as seen in the previous 
paragraphs, EDE manifests itself not only in the shift of the 
peaks, but also in the shape of the peaks and in the low-Z ISW 
effect. In this example, the height of the first peak is differ- 
ent for the EDE model, as is the low I behaviour, rendering it 
distinct from the flat ACDM model. Thus, although we expect 
some degeneracy between the dark energy parameters and the 
curvature, this degeneracy is not very strong, since both the po- 
sition and the height of the first peak are strongly constrained 
by current CMB data. 

3.2. Constraints from Observations 

We first study the results using only the CMB data. 
The primary parameters to be varied are the standard CMB 
parameters- fl c h 2 ,flt,h 2 ,8,T,n s ,A s , and the equation of state 
parameters wo,w m ,a,, A,. Since we wish to study EDE 
models, we restrict the equation of state at early times to 
w m < 0.1. This can give rise to a dark energy density of up 
to Hde(z) < 0.03f2 m (z) at early times. We assume a flat uni- 
verse, i.e. O re = 0, and consider the full dark energy perturba- 
tions. The secondary parameters that we deduce from the anal- 
ysis are Qom,Ho, fW/^mCzcm))^ The first column of ta- 
ble [TJ shows the mean and 2a boundaries for the primary and 
secondary parameters. We see that the EDE parameters are 
constrained at wq < -0.61, a, < 0.44 (which means z r > 1.2), 
A, < 0.37. The constraints on the scale and width of transi- 
tion are reasonable, however, the constraint on the equation of 
state today, Wq, is too broad. We note however, that this result 
is obtained by using CMB data alone, using other data would 
reduce degeneracies with the other parameters. For instance, 
SNe Type la data would affect the equation of state today more 
strongly. Also, the Hubble parameter for which wq ~ -0.6 is 
allowed is Hq ~ 60 km/s/Mp c, much lower th an the currently 
accepted measurement for it (Rie ss et al. 20091) . Therefore, we 
expect that the addition of other observations to the analysis 
should improve the constraints on the EDE parameters signifi- 
cantly. It is interesting to note also that the 2a lower bound on 
erg for this analysis is as low as a& > 0.49, whereas the ACDM 
fit to the WMAP5 data has a a & > 0.72. 

We now redo the analysis adding other datasets to see how 
the constraints improve. Three distinct cases are considered- 
(a) full dark energy perturbations are taken into account, fl K = 
0; (b) dark energy perturbations are considered negligible, Ct K = 
0; and (c) full dark energy perturbations are considered, and 
the constraint on the flatness of the universe is lifted (i.e. £l K 
is a free parameter). The results are shown in the second, 
third and fourth columns of table [TJ For the fully perturbed, 
flat case, when all the data is considered, the EDE parame- 
ters are constrained to wq < -0.89, a, < 0.19 (i.e. z, > 4.2), 
A, < 0.21. The addition of other datasets clearly enhances 
the constraints on the EDE model. This is because the other 
parameters which could be degenerate with the EDE param- 
eters, such as the Hubble parameter, are well-constrained by 
other observations. We note that the constraint on the equation 
of state today, Wo, is stronger than that would be obtained us- 
ing the background data alone (e.g., for Type la SNe, we find 
wo <, —0.75 for constant equatio n of state, when systematics 
are included dHicken et al. 2009)). In figure [4] we show the 
two-dimensional 68% and 95% confidence levels, as well as 
the marginalized one-dimensional distributions for the EDE pa- 
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rameters of interest, wo,a t ,A t , and the matter density Slo m and 
the Hubble parameter Hq, which are expected to be degenerate 
with the EDE parameters. We see that all three EDE parameters 
are now strongly constrained, and the non-EDE parameters are 
close to the values expected in the ACDM model. The evolution 
of the equation of state of dark energy with redshift is shown in 
figure|5] We see that at low redshifts (z _< 2), the 2a confidence 
level for w(z) is quite close to ACDM. Thus current observa- 
tions already constrain the evolution of the equation of state for 
EDE models to ACDM-like behaviour at present and in the near 
past. Studying the background expansion data (which is usu- 
ally below redshift of two) will therefore not able to distinguish 
these EDE models from ACDM with any success even if there 
is adequate amounts of dark energy at early times. In order to 
distinguish these EDE models (currently accepted by the data) 
from ACDM, we need to look at the perturbative data. Thus we 
may conclude that even if there is significant amount of dark 
energy in the universe at early times, this has to reduce to dark 
energy very close to ACDM at present times, that this transition 
cannot take place too late (around redshift of four) and that the 
transition needs to be sharp (A, < 0.2). However, we should 
note that even with these constraints, erg is still significantly dif- 
ferent from the typical ACDM value, with the 2a lower bound 
being at erg ^> 0.6. This is because, as discussed earlier in sec- 
tion l3.ll EDE has a strong effect on the matter power spectrum, 
leading to a much lower erg than that in the cosmological con- 
stant model. This means that studying data which utilizes the 
matter power spectrum (such as the halo mass functions) even 
at low redshifts may allow us to discriminate between EDE and 
ACDM models. 

We next look at the case where dark energy perturbations are 
neglected, for the full dataset. We find that, although the results 
are similar for many of the parameters, they can be rather dif- 
ferent for the EDE parameters. As seen in the third column of 
table Q] the EDE parameters are constrained to wq < -0.8, a, < 
0.33 (i.e. Zt > 2), A, < 0.31, thus, neglecting the perturbations 
for an EDE model would result in rather broader constraints 
on its parameters. The value of erg allowed at 2a is also much 
closer to the ACDM value, with erg > 0.72. Neglecting the dark 
energy perturbations in a dynamic dark energy scenario may 
therefore produce results very different from the true results 
when full dark energy perturbations are considered. 

If the flatness of the universe constraint is removed, taking 
dark energy perturbations into account, the EDE parameters 
are mildly degenerate with the curvature of the universe 51 K . 
As seen in the last column of table Q] the EDE parameters in 
this case are constrained to w$ < -0.77, a, < 0.35 (i.e. z, 2), 
A, < 0.35, while the curvature of the universe is still rather 
tightly constrained to -0.014 < Q, K < 0.031. Thus, relaxing the 
flatness constraint leads to a weakening of the constraints on the 
parameters of the EDE models, but can still lead to reasonable 
constraints on the EDE parameters. 

Previous works that have studied EDE with perturbations 
have constrained th e amount of early dark energy using current 
observations, e.g., dXia & Viel 20 09) obtained VI E de < 1.4 x 
10~ 3 . However, as explained in section [2~2l this study allowed 
for models of dark energy that have negligible amounts of dark 
energy at early times. This led to a weakening of the constraints 



on the EDE transition parameters, with the parameter w m attain- 
ing peaks both at ~ -1 and ~ 0, and the equation of state today 
being close to ACDM. Thus no strong constraint could be put 
on the evolution of the equation if state. In this work we have 
attempted to address the question of how to put constraints on 
the transition of early time to late time dark energy if the uni- 
verse contains a certain amount of early dark energy. If we 
constrain early dark energy to w m >, -0.1, we exclude models 
which do not have EDE behaviour, and thus are able to put rea- 
sonable constraints on the transition parameters, which give us 
an insight into the evolution of dark energy for these models. It 
is difficult to rule out the presence of EDE altogether, due to the 
dearth of data at very high redshifts, but with this study we are 
able to put constraints on when the universe could have tran- 
sited from such early time dark energy to late time, ACDM-like 
behaviour. We find that late time behaviour of the equation of 
state of these models must be close to ACDM below redshift 
of few. However, since the erg of these models is rather differ- 
ent from ACDM, they may be distinguished from ACDM using 
data such as the halo mass function, or the Sunyaev-Zeldovich 
power spectrum, even at the lower redshifts. In addition, this 
work also studies the degeneracy between the curvature of the 
universe and the dark energy parameters. 

4. CONCLUSIONS 

In this work, we have studied early dark energy models us- 
ing current observations. We find that, if a sizeable amount of 
dark energy exists in early times (Hde(zcmb) — 0.03il m (zcMB)), 
we may put tight constraints on the transition of this dark en- 
ergy to its present day value, and that the present day value of 
the dark energy equation of state must be close to the ACDM 
value. If the dark energy perturbations are correctly accounted 
for, then the current dark energy equation of state is con- 
strained to Wo < -0.89, while the transition from early dark 
energy must occur at redshifts of Zt > 4.2, with a narrow transi- 
tion width of A, < 0.21. Incorrectly assuming that dark en- 
ergy perturbations are negligible leads to a different result- 
wq < — 0.8,z ( > 2, A, > 0.31, thus showing that it is vital to 
include the dark energy perturbations self-consistently in any 
analysis that uses perturbative data such as CMB or the matter 
power spectrum. Leaving ft K to be a free parameters leads to 
a weakening of the constraints on the dark energy parameters, 
with w < -0.77, a, < 0.35 (i.e. z, > 2), A, < 0.35 for -0.014 < 
fi K < 0.03 1 . We note that, for the flat universe in which dark en- 
ergy perturbations are considered, the value of erg is much lower 
than that in correspo nding ACDM models. As will b e shown in 
a companion paper (lAlam et al. in preparation 20101) . this may 
lead to interesting constraints from future large scale structure 
data such as halo mass functions, as also from the Sunyaev- 
Zeldovich power spectrum. 
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Fig. 1 . — Growth of relative dark energy density Pde/ Pm (panels (a), (d)), dark energy density contrast Sde & velocity perturbation \>de (panels 
(b), (e)), and cold dark matter density contrast 5 m (panels (c), (f)) with the scale factor a for two dark energy models : DE with wde = -0.9, 
and EDE with wo = -1.0, w,„ = -0.1, a, = 0.3, A, = 0.2 respectively. The red line in panels (a) and (d) represents the DE and EDE models 
respectively, the black line in panels (a), (c), (d), (f) shows the ACDM model for comparison, the dark energy model without perturbation is 
shown in green in panels (c) and (f), while the case with perturbation is shown in red in these panels. The solid line in panels (b) and (e) represents 
the dark energy density contrast, while the dashed line shows the velocity perturbation. The filled circle in panels (a), (d) represent matter-dark 
energy equality for the DE and EDE models respectively, while the cross represents matter-dark energy equality for ACDM. 
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Fig. 2. — Scalar C/'s and matter power spectrum for the DE model with wde = -0.9 (panels (a) and (b)), and the EDE model with wo = -1.0, w m = 
-0. 1 , a, = 0.3, A,- = 0.2 (panels (c) and (d)). The black line in each panel represents the corresponding ACDM model, the green line represents 
the dark energy model with no dark energy perturbations, while the red line represents the case with dark energy perturbations taken into account. 
The insets in the panels (a) and (c) show the shift in the position of the first peak for the dark energy model considered. Panels (e) and (f) show the 
difference between the perturbed and unperturbed cases for both the DE (solid line) and the EDE (dashed line) models for the two observables. 
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Fig. 3. — Scalar Q's for ACDM and EDE models with curvature. The black line in each panel represents the flat ACDM model, the green line 
represents ACDM with fi K = 0.06, and the red line represents an EDE model with wo = -0.65, w m = -0.1, at = 0.2, At = 0.1, with fl K = 0.06. The 
inset shows the shift in the position of the first peak for the ACDM model with curvature. For the EDE model with curvature, the shift in the first 
peak is compensated, however, the height of the first peak as well as the low-/ behaviour differs from the flat ACDM case. 
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Fig. 4.— Constraints from CMB (WMAP5, CBI, VSA, BOOMERANG, ACBAR) and other datasets (SNe Type la Constitution, LRGDR7, 
SHOES) on EDE parameters wo, a,, A,, as well as fiom,#o, showing marginalized one-dimensional distributions and two-dimensional 68% and 
95% limits. Full dark energy perturbations are taken into account and the curvature of the universe is fixed at Q K = 0. 
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Fig. 5.— 2a confidence levels from CMB (WMAP5, CBI, VSA, BOOMERANG, ACBAR) and other datasets (SNe Type la Constitution, 
LRGDR7, SHOES) on the equation of state of dark energy, for the EDE models considered. Full dark energy perturbations are taken into account 
and the curvature of the universe is fixed at Q K = 0. 



